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ABSTRACT 


Ifc 


Density  functional  formalism  of  Hohenberg  and  Kohn  (HK)  is  generalized 

for  the  case  of  a multicomponent  plasma.  Using  the  self-consistent 

K-  hn-Sham  (KS)  equations  for  electrons  and  holes  and  local 

density  approximation  for  exchange-correlation  potential , 

we  investigate  the  surface  characteristics  of  electron- 

hole  liquid  (EHL)  in  six  configurations  of  Ge  and  SI . We  denote  these 

configurations  by  X(v  ;v.), where  X is  either  Ge  or  Si,  and  v and  v are  the 

c n 7 e n 

number  of  occupied  electron  and  hole  bands* respectively.  In  normal  Ge,  i.e., 

—4  2 

Ge(4;2),  the  value  of  surface  tension,  o,  is  found  to  be  3.7  x 10  erg/cm  . 

2 

When  Ge  is  subject  to  a uniform  stress  of  about  3.5kg/mm  along  <111>  direc- 

-4  2 

tion,  i.e.,  in  Ge(l;2),  o is  calculated  to  be  1.0  x 10  erg/cm  . Under  a 
very  large  <111>  uniaxial  stress  on  Ge,  i.e.,  Ge(l;l),  o is  found  to  be  a fac- 
tor of  twenty  smaller  than  in  Ge(4;2). 

Charge  on  electron-hole  drop  (EHD)  is  also  studied  in  the  above  mentioned 
systems.  In  accordance  with  the  experiment  of  Pokrovsky  and  Svlstunova,  we 
find  that  EHD  is  negative  in  Ge(4;2)  and  positive  in  Ge(l;2).  It  is  predicted 
that  the  drop  will  sustain  a negative  charge  in  Ge(l;l). 

Calculations  for  surface  tension  and  charge  on  EHD  are  also  reported  in 

three  configurations  of  silicon.  The  value  of  a in  unstressed  Si,  denoted  by 

4 2 

Si(6;2),  is  obtained  to  be  87.4  x 10  erg/cm  . Application  of  an  intermedi- 
ate stress  along  <100>  direction  leads  to  the  configuration  Si(2;2).  The 

_4  2 

value  of  a in  Si(2;2)  is  found  to  be  32.8  x 10  erg/cm  . In  the  presence  of 
a large  <100>  stress,  i.e-,  in  Si(2;l),  surface  tension  is  a factor  of  eight 
smaller  than  in  Sl(6;2).  Calculation  of  charge  reveals  that  the  EHD  is  nega- 
tive in  both  Si(6;2)  and  Si(2;l).  Within  the  limits  of  accuracy  of  our  calcu- 
lation we  find  the  drop  is  almost  neutral  in  Si(2;2). 
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INTRODUCTION 


In  indirect  band  gap  semiconductors  (e.g.,  Ce  or  Si)  excitons  are  formed 
under  weak  external  excitation.  When  the  density  of  excitons  becomes  high, 
their  individuality  is  lost.  Further,  in  the  luminescence  spectra  a new  line 
appears  on  the  low  energy  side  of  the  exciton  line.  Detailed  experimental  in- 
vestigations have  consistently  pointed  out  that  the  new  luminescence  line  ori- 


ginates from  a recombination  of  electrons  and  holes  condensed  in  high  density 


droplets^  Such  a condensed  state  is  now  known  as  electron-hole  liquid 


(F.HL) . Experimental  investigation  of  its  properties  is  facilitated  by  the 

,-6 


fact  that  the  life  time  of  particles  is  quite  long  (~  10  sec)  and  therefore 
the  condensed  phase  is,  to  a very  good  approximation,  in  a state  of  thermo- 
dynamic equilibrium.  It  was  Keldysh  who  first  conceived  that  the  EHL  is  lik- 


ened to  a plasma  of  electrons  and  holes  stabilized  by  Coulomb  interaction  of 


.(17) 


the  constituents  . His  conception  of  EHL  was  subsequently  confirmed  by 


an  abundance  of  experimental^  and  theoretical^^  ^ observations. 


Theoretical  appeal  in  electron-hole  liquid  is  mainly  due  to  the  absence 
of  complications  which  can  arise  from  the  lattice  effects.  In  view  of  the 

• o 

fact  that  the  excltonic  Bohr  radius  in  Ce  (-  177  A)  or  Si('49  A)  is  large 
compared  to  the  lattice  constant,  and  that  the  fraction  of  electrons  excited 


into  conduction  band  is  small,  the  interaction  between  electrons  and  holes  is 

(27) 


bare  Coulomb  screened  by  static  dielectric  constant  of  the  material  . Thus, 
EHL  can  be  regarded  as  a collection  of  electrons  and  holes,  characterized  by 
proper  band  masses  and  interacting  via  statically  screened  Coulomb  interaction. 
The  fact  that  the  system  at  hand  can  be  accurately  described  as  a quantum  plas- 
ma of  electrons  and  holes  makes  it  an  ideal  candidate  for  testing  different 
many-body  approximations.  Theoretical  efforts  to  date  have  been  successful  in 
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( 20-26} 

explaining  the  ground  state  properties  of  EHL'  . Calculations  of  bind- 

ing energy  and  equilibrium  density  in  Ge  and  Si  agree  very  well  with  experi- 
ment'  . Detailed  investigations  have  also  been  made  of  thermodynamic  quan- 
tities like,  gas-liquid  transition  temperature,  the  compressibility,  the 

(16  26  28-33) 

temperature  dependence  of  density  and  of  chemical  potential,  etc.  * ’ 

Theoretical  results  for  these  quantities  are  also  in  good  agreement  with  ex- 

( 26} 

periments,  thereby  establishing  the  validity  of  plasma  model' 

Surface  structure  of  EHL  has  not  hitherto  been  examined  as  thoroughly  as 
its  bulk  properties.  The  reasons  for  insufficient  knowledge  of  EHL  surface 
are  essentially  two  fold.  First,  the  lack  of  translational  invariance  makes 
any  surface  problem  more  difficult  than  the  bulk.  Second,  proper  understand- 
ing of  the  bulk  properties,  which  must  precede  that  of  the  surface,  has  come 
about  only  recently.  Now  that  the  bulk  properties  of  EHL  are  well  understood, 
there  is  a basis  for  a proper  investigation  of  its  surface  properties. 

In  this  paper  we  shall  examine  the  surface  characteristics,  such  as  sur- 
face tension,  dipole  layer  and  charge  on  electron-hole  droplet  (EHD)  in  ger- 
manium and  silicon  under  uniform,  uniaxial  stress  along  <111>  and  <100>  direc- 
tions, respectively.  For  the  sake  of  convenience  we  shall  designate  these 
(34) 

systems  by  X(ve;v^)  where  X is  either  Ge  or  Si,  and  and  vh  are, respec- 
tively, the  number  of  conduction  and  valence  bands.  In  germanium  under  zero 
stress,  i.e.,  Ge(4;2),  there  are  four  equivalent  conduction  bands  along  <111> 
direction  and  two  hole  bands  degenerate  at  the  T point.  Away  from  the  center 
of  Brillouin  zone,  the  hole  bands  branch  out  into  a light  and  heavy  hole  bands. 
The  structure  of  valence  bands  is  the  sam*  in  unstrained  silicon  and  germanium. 
The  difference  arises  in  the  number  and  location  of  conduction  valleys.  There 
are  six  conduction  bands  along  <100>  direction  in  Si.  Thus,  the  notation  for 
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unstrained  silicon  is  Si(6;2; 

In  the  presence  of  a <111>  stress  on  Ge, one  of  the  conduction  minima  moves 

toward  the  valence  bands  while  the  remaining  conduction  bands  move  away  from 

2 

it.  At  a stress  of  about  3.5kg/tm  all  the  electrons  of  the  condensed  phase  re- 
side in  the  lowest  conduction  valley*  because  the  Fermi  energy  of  electrons  co- 
incides with  the  absolute  minima  of  the  other  three  conduction  bands.  Except 
for  a slight  decoupling  at  the  T point,  the  structure  of  valence  bands  remain 
as  in  Ge(4;2).  Such  a configuration  is  called  Ge(l;2).  When  silicon  is  subject 
to  a uniform  stress  along  <100>  direction,  two  of  the  six  conduction  valleys 
move  toward  the  valence  bands  while  the  remaining  four  move  away  from  it.  At 


a stress  of  about  10.5  kg/nrni  , 


only  the  two  lowest  conduction  bands 


are  populated.  The  changes  that  occur  in  the  structure  of  valence  bands  are 
similar  to  those  in  Ge(l;2).  Following  our  notation  we  denote  this  configu- 
ration by  Si(2  ;2) . 

Application  of  a large  <111>  uniaxial  stress  in  Ge  and  <100>  in  Si  not 

only  removes  the  valence  band  degeneracy  at  the  T point,  but  also  modifies  the 

(36) 

structure  of  valence  bands.  In  both  systems,  the  holes  occupy  a single,  highly 
anisotropic  band.  In  Ge  under  a large  <111>  uniaxial  stress  the  electrons  re- 
side in  a single  conduction  band,  whereas  in  Si  under  a large  <100>  stress  they 
occupy  two  conduction  bands.  We  designate  these  systems  by  Ge(l;l)  and  Si(2;l). 
The  band  masses,  dielectric  constant,  and  the  values  of  excltonic  rydberg  used 
in  the  present  work  are  listed  in  Tables  I and  II. 

Owing  to  different  masses  and/or  number  of  bands,  electrons  and  holes  in 
EHL  possess  different  chemical  potentials.  Now,  the  bulk  chemical  potential  for 
any  component  consists  of  kinetic,  exchange  and  correlation  contributions.  Ex- 
plicit calculations  for  the  ground  state  energy  reveal  that  electrons  and  holes 


I 


contribute  almost  equally  to  exchange-correlation  energy 


(21-26) 


In  that  event. 


the  difference  in  bulk  chemical  potentials  of  electrons  and  holes  arises  mainly 
— T —X  -T  -t 

from  (u  - p.),  where  p and  p,  are,  respectively,  the  Fermi  energies  of  elec- 
eh  eh 

trons  and  holes.  For  a system  with  vg  conduction  bands,  and  a light  and  heavy 
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hole  bands,  v and  p.  are  given  by 
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where  nQ  is  the  equilibrium  density  of  e-h  pairs,  m^e  the  density  of  states 
mass  for  an  electron,  and  m^  and  are,  respectively,  the  light  and  heavy 
hole  masses. 

-T  — T 

Using  the  band  masses,  given  in  Table  I,  for  Ge(4; 2)  we  find  that  jj^/ ■ 
2/3 

1.62/v  - 0.64,  which  implies  that  the  holes  in  EHL  are  less  tightly  bound 

than  the  electrons.  Consequently,  at  a low  but  finite  temperature  the  holes 
will  have  excess  thermionic  emission  and , therefore , the  EHD  will  sustain  a 
negative  charge.  The  holes  will  continue  to  evaporate  until,  at  equilibrium, 
the  work  functions  for  electrons  and  holes,  and  therefore  their  rates  of  emis- 
' ' . Similarly,  in  Si(6;2)  the  ratio  p*/p^  ■ 0.55,  imply- 


sion  become  equal 
ing  that  the  EHD  will  be  negatively  charged. 

Application  of  these  arguments  leads  to  interesting  consequences  in 
strained  Ge  and  Si.  It  is  easy  to  infer  that  in  Ge(l;2)  and  Si(2;2)  the  elec 

1.62  in 


-T  -T 

trons  will  tend  to  evaporate  more  than  the  holes  (the  ratio  p /p, 

e n 


Ge(l;2)  and  1.14  in  Si(2;2))  and  therefore  the  droplet  will  acquire  a net  posi- 

-T  -T  2/3 

tive  charge.  In  the  case  of  Ge(l;l)  and  Si (2; 1)  the  ratio  y^/p^  ■ mdj/ndeVe 
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With  the  band  masses  given  in  Table  II,  “ 0-40  in  Ge(l;l)  and  0.46  in 

Si(2  ;l) , which  implies  that  the  electrons  will  evaporate  less  readily  than  the 

holes.  In  other  words,  the  EHD  will  be  negatively  charged  in  Ge(l;l)  and  Si(2;l). 

Measurements  of  Pokrovsky  and  Svistunova  have  confirmed  that  the  charge  on  EHD 

(38) 

is  indeed  negative  in  Ge(4;2)  and  positive  in  Ge(l;2)v  ' . To  date,  no  experi- 
mental measurement  of  droplet  charge  has  been  reported  in  Ge(l;l),  Si(6;2), 

Si(2;2)  and  Si(2;l). 

Coming  to  the  question  of  surface  tension,  experimental  studies  of  coexist- 
ence curve  have  revealed  that  the  exciton  gas  becomes  supersaturated  before  the 

(39) 

EHD  begins  to  nucleate  . This  is  a signature  of  existence  of  surface  tension. 

Measurements  of  EHL  surface  tension,  a,  have  been  carried  out  only  in  Ge(4;2). 

-4  2 (39) 

Westervelt  et  al.  have  reported  a value  of  2.9*10  erg/cm  . Bagaev  and  co- 

-4  2 (40) 

workers  have  estimated  o to  be  1.6*10  erg/cm  . Recent  experiment  of  Eti- 

enne et  al.  reveals  a value  of  3 * 10  ^erg/cm^ while  Staehli's  estimate  is 
3.8  * 10  ^erg/cm2  ^2\  To  our  knowledge,  there  does  not  exist  any  experimental 
measurement  of  surface  tension  in  Ge(l;2),  Ge(l;l),  Si(6;2),  Si(2;2),  and  Si(2;l). 
Without  conclusive  experiments  on  surface  tension  and  binding  energy  it  is  in- 
conceivable that  any  progress  can  be  made  in  understanding  the  coexistence  curve 
or  the  kinetics  of  EHD  formation. 

In  view  of  the  fact  that  EHL  is  free  from  "parasitic"  effects  of  ions,  theo- 
retical investigation  of  its  surface  properties  can  be  appropriately  handled  by 

(43-45) 

means  of  density  functional  formalism  of  Hohenberg-Kohn-Sham  (HKS)  . With- 

in this  formalism  there  are  essentially  two  ways  of  approaching  the  surface  prob- 
lem. The  easier  of  these  two  is  widely  known  as  variational  procedure.  The  basic 
philosophy  of  variational  method  is  to  approximate  the  ground  state  energy  and 
use  exponential  density  profiles  for  electrons  and  holes.  Each  density  is  char- 
acterized by  a variational  parameter.  In  this  way  the  approximate  ground  state 


-6- 


energy  becomes  a function  of  variational  parameters.  After  minimizing  total 
energy  with  respect  to  these  parameters  one  obtains  the  minimized  total  en- 
ergy. Subtracting  the  bulk  contribution  from  the  total,  one  gets  the  surface 
energy.  Most  authors  have  retained  local  density  contributions  from  kinetic 
and  exchange-correlation  energies,  and  first  gradient  correction  to  kinetic 

energy^^’^  The  value  of  surface  tension  that  they  obtain  in  Ge(4;2) 

-4  2 

is  around  1*10  erg/cm  . On  including  the  first  gradient  correction  to  ex- 

change-correlation  energy  Vashishta,  Kalia  and  Singwi  find  that  the  value  of 

surface  tension  changes  significantly (50,51),  Their  estimate  of  o in  Ge(4;2) 
-4  2 

is  3.5*10  erg/cm  . Evidently,  it  is  in  reasonable  agreement  with  the  meas- 

(39)  (41)  (42) 

urement  of  Westervelt  et  al.  , Etienne  et  al.  , and  cf  Staehli  . An 

important  outcome  of  their  variational  calculation  as  well  as  that  of  Reinecke 
(49) 

et  al.'s  is  that  the  EHD  surface  tension  in  Ge(4;2)  is  an  order  of  magni- 
tude (a  factor  of  sixteen)  larger  than  in  Ge(l;l)/50^ 

/ / Q \ 

Variational  calculations  that  include  effects  of  valence  band  coupling'' 
on  gradient  correction  to  kinetic  energy  of  holes,  as  well  as  exchange-corre- 
lation gradient  correction , yield  a positive  charge  on  EHD  in  Ge(4;2)  and 
negative  charge  in  Ge(l!2).  These  results  are  in  disagreement  with  the  above 

mentioned  thermodynamic  arguments  and  the  experimental  results  of  Pokrovsky 
(38) 

and  Svistunovav  ' . Further,  contrary  to  what  one  expects  from  energetic  con- 
siderations, the  variational  procedure  yields  a positive  charge  on  EHD  in 
Ge(l;l)^^.  Such  an  unsatisfactory  feature  of  variational  method  arises 
from  the  approximate  treatment  of  kinetic  energy  and  the  absence  of  Friedel 
oscillations. 

Considerable  improvement  can  be  made  over  variational  procedures,  provided 

(44  45) 

one  resorts  to  the  self-consistent  scheme  of  Kohn  and  Sham  ’ ; the  work 

involved  in  it  is  stupendous.  The  merit  of  self-consistent  procedure  lies 
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in  the  exact  treatment  of  kinetic  energy  which  also  brings  in  part  of  the  ef- 
fect of  Friedel  oscillations.  We  shall  describe  in  this  paper  a self-consist- 
ent calculation  for  the  surface  properties  of  EHL  in  the  aforementioned  sys- 

(38) 

terns.  In  accordance  with  the  observation  of  Pokrovsky  and  Svlstunovav  , the 

self-consistent  calculation  yields  a negative  charge  on  EHD  in  Ge(4;2),  and  a 

(52) 

positive  charge  in  Ge(l;2)  . On  the  basis  of  our  calculation  we  predict 

that  the  EHD  will  be  negatively  charged  in  Ge(l;l),  Si(6;2),  and  Sl(2;l).  This 

is  in  agreement  with  the  aforementioned  thermodynamic  considerations.  Within  the 

limit  of  accuracy  of  the  calculation,  the  EHD  is  found  neutral  in  Si(2;2). 

-4  2 

As  regards  surface  tension,  we  obtain  a value  of  3.7  * 10  erg/cm  in 

Ge(4;2)  which  is  in  reasonable  accord  with  the  experimental  estimates  of  West- 

09)  (41)  (42) 

ervelt  et  al.  , Etienne  et  al.  , and  Staehli  . Surface  tension  in 

GeO  »1)  is  found  to  be  a factor  of  twenty  smaller  than  in  Ge(4;2).  The  pre- 

-4  2 

sent  calculation  yields  a value  of  87.4  * 10  erg /cm  for  surface  tension  in 

Si(6;2).  As  in  the  case  of  Ge  under  stress,  we  find  that  the  o in  Si(2;l)  is 

much  smaller  (=  1/8)  than  in  Sl(6;2).  In  Ge(l;2)  and  Si(2;2)  the  values 

-4  -4,2 

of  a are  calculated  to  be  1.0  * 10  and  32.8  * 10  erg/cm  , respectively. 

In  Section  I we  give  a multicomponent  generalization  of  HK  theory  and  de- 
rive rigorous  expressions  for  the  chemical  potential  of  electrons  and  holes. 

It  is  followed  by  a description  of  self-consistent  Kohn-Sham  equations  in 
Section  II.  Section  III  deals  with  the  surface  tension  of  EHD  in  the  six  sys- 
tems and  Section  IV  is  devoted  to  a discussion  of  EHD  charge.  Finally,  we 
mention  in  the  conclusion  vhat  improvements  can  be  made  on  the  present  calcu- 
lation. 

I.  MULTICOMPONENT  GENERALIZATION  OF  HK  THEORY 

Consider  an  M-component  system  in  the  presence  of  external  potentials 
(V*(r)},  each  of  which  couples  with  only  a particular  component.  The  Hamil- 
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tonlan  of  such  an  M component  system  can  be  expressed  as 


M 

H“i£i 


d—  <-  h v2) 


ft(r) 


M 

i£i 


| dl  ** (r)  V*(r)  ¥t(r) 


(2) 


+ 


M M 

1/2  ih  ih 


where  ¥^(r)  is  a field  operator  for  the  ic^  component,  ¥^(r)  its  Hermitian 
conjugate,  and  v„(r)  is  the  statically  screened  Coulomb  potential 

energy.  The  expectation  value  of  H taken  with  respect  to  the  ground  state 
wavefunction  |¥>  yields  the  ground  state  energy,  E: 


E - T + V + e J dr  V*(r)  n^r)  , (3) 

where  » + 1 for  holes  and  -1  for  electrons.  In  Eq.  (3),  T and  V are  the 
kinetic  and  potential  energies,  respectively,  and  n^(r)  is  the  density  of  the 
i.th  component. 

(43) 

Following  the  proof  given  by  HK  for  a single  component  system  , we  shall 
show  by  reductio  ad  absurdum  that  the  ground  state  energy  of  a multicomponent 
system  is  a unique  functional  of  the  density  of  each  component.  Let  us  assume 
that  there  exists  two  sets  of  external  potentials,  (V*a(r)}  and  {V*g(r)},  and 
that  corresponding  to  them  there  are  two  wavefunctlons  |v>  and  |¥’>  whif'h  yield 
the  same  set  of  densities  (n^(r)}.  We  shall  further  assume  that  the  ground 
state  is  nondegenerate  and  that  the  ground  state  energy  is  minimum  with  respect 
to  variations  in  wavefunction.  In  that  case. 


<«#» |h|4»)  < - (*'|hV)  + (♦'Ih  “ (4) 

or 

E < E'  + ifi  e ?i  | dr  n^r)  [V^(r)  - V*g(r)]  (5) 

Interchanging  the  primed  and  unprimed  quantities  and  adding  the  resulting 
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equation  to  Eq.  (5),  we  obtain 


I 


E + E'  < E + E'  , 


(6) 


which  is  absurd.  Therefore,  |f>  and  cannot  be  different  and  yet  yield 

the  same  set  of  densities  (n^(r)}.  In  other  words,  the  external  potentials 

{ (r) } and,  therefore,  the  ground  state  energy  are  uniquely  determined  by 
e — 

If  the  number  of  particles  is  fixed,  E must  be  stationary  with  respect 
to  variations  in  any  n^(r).  A change  in  n^(r)  from  its  correct  value  shall 
correspond  to  a change  in  the  wavefunction  from  the  ground  state  function 
and  according  to  variational  principle  such  a state  is  of  higher  energy. 

Insofar  as  the  potential  energy  has  a long  range  Coulomb  part,  it  is 
hardly  convenient  to  express  the  ground  state  energy  as  in  Eq.  (3),  i.e., 
a sum  of  contributions  arising  from  external  perturbation,  potential  and 
interacting  kinetic  energies.  Instead,  it  is  much  more  appropriate  to  sep- 
arate the  classical  Coulomb  energy  and  write  the  ground  state  energy  of  a 
multicomponent  plasma  as 


z m M f Ml' 

E t{ni»  - b ih  ih  j d*’  d^"  ir'  --Vr 


(7) 


+ 6 ill  «i  | Ve(i’>  ni(£’>  + ih  Ts[ni]  + Exc[{ni}]  * 


where  the  summations  are  taken  over  different  components  of  EHL.  The  first 
term  in  Eq.  (7)  constitutes  the  classical  Coulomb  energy,  the  third  term  is 
the  sum  of  noninteracting  kinetic  energies  of  each  component  and  the  last 
term  is  the  exchange-correlation  contribution.  A distinct  advantage  of  writ- 
ing the  ground  state  energy  in  a form  such  as  Eq.  (7)  lies  in  the  ease  with 
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which  different  terms  in  the  equation  can  be  expressed  in  terms  of  quantities 
related  to  the  homogeneous  system. 

Let  6nt  denote  an  arbitrary  variation  in  the  density  of  component. 
Owing  to  the  stationary  property  of  E, 

6JtE[{n±>]  - 0 . (8) 

Since,  the  number  of  particles  of  each  component  is  conserved. 


djr  6n^(r)  ■ 0 


(9) 


Conditions  expressed  in  Eqs.  (8)  and  (9)  are  basic  to  density  functional 
formalism.  Combining  them  by  means  of  an  undetermined  Lagrange  multiplier, 
we  obtain 


^Eftn^]  - X di  ^(oj  - 0 


Using  Eq*  (7)  in  Eq.  (10)  we  obtain  the  fallowing  expression  for 

2,  # l \ •> i 


Ut  - e e 


I e%  M f n^r')  61%  In]  «E  [{n » 

ive(r>  + -T1  Jll  I ^ \^T\  + -l^T  + 6n'(r)— 


(10) 


(ID 


Henceforth,  we  shall  suppress  the  external  potential. 

Let  £^(r)  represent  the  electrostatic  energy  of  component.  With 
<Kr)  given  by 


e2  M 
~ j-1  5j  J 


n.(r') 
dr'  — * 


- I r - r' 


(12) 


the  expression  for  reads. 


wi  " + «nt(r) 


6W  , 6Exct{ni}l 


6njl(r) 


(13) 


-4 

In  view  of  the  fact  that  the  radius  of  EHD  (*  5 x 10  cm)  is  much  larger 


than  interparticle  distance  in  the  system,  it  is  reasonable  to  assume  that 


EHL  is  a semi-infinite  system,  occupying  half  the  space  z < 0.  The  space  z > 0 
is  treated  as  vacuum.  Under  these  conditions,  densities  of  EHL  components  vary 
only  in  the  z direction,  thus  making  4 a function  only  of  z.  Transforming  Eq.  (12) 
into  a Poisson  equation,  we  obtain 


d24(z)  4ve2  £ , / * 


(14) 


The  appropriate  boundary  conditions  for  4(z)  are  taken  to  be, 


*(z)  z - + 


constant 


and 


dz 


(15a) 


(15b) 


'z  -*■  + - 

Then,  the  solution  of  Poisson  equation  reads. 


4ire2  M C 

4(*)  - 4(-)  + — Cj  j dz ' (z  - z’)nj(z')  (16) 


Following  the  conventional  definition  of  dipole  layer,  namely 


64  - 4(-)  - 4 (-  -)  , (17) 

and  making  use  of  the  fact  that  the  total  system  is  charge  neutral,  we 

obtain  2 /•» 

j£i  j « 4*'  *'  nj<*>  U«> 


Writing  Eq.  (18)  in  terms  of  electron  and  hole  densities,  ne(z)  and  n(,(z),  we 

get  4ve2  r 

A4  - - dz'  z'  [nt(z')  - nh(z')]  (19) 

It  is  consequential  to  the  discussion  of  charge  to  express  the  chemical 
potential  in  terms  of  dipole  layer.  For  this  purpose,  v*  average  Eq.  (13) 
over  the  volume  of  the  drop  and  obtain 


r I(  , , , SE,c»"l>' 

~w. — 


(20) 
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where  is  the  density  of  tch  component  in  the  bulk.  The  last  two  terms  in 
Eq.  (20)  are,  respectively,  the  kinetic  and  exchange-correlation  contributions 
to  bulk  chemical  potential  of  l1-*1  component.  The  expression  for  kinetic  con- 
tribution Is  simply  the  Fermi  energy  for  that  component.  The  exchange-corre- 
lation part  of  chemical  potential  is  obtained  from  the  ground  state  energy 

(23  26) 

calculation  for  a uniform  systemv  ’ ' . 

From  now  onwards  we  shall  take  $(<»)  * 0 as  our  reference  level  of  energy. 
In  that  case,  Eq.  (20)  becomes 


«T4Ca.3  6E  {{n.}] 
s t + xc  i 


<5Ha 


- V* 


(21) 


The  discussion  of  charge  that  we  presented  in  the  introduction  can  now  be 
made  completely  general  provided  one  compares  the  chemical  potential,  of 
electrons  with  that  of  holes,  y^.  Since  the  bulk  chemical  potential  is  ob- 
tained from  the  ground  state  energy  of  a uniform  EHL,  the  quantity  of  main 
interest  in  a surface  calculation  is  the  dipole  layer.  The  dipole  layer 
depends  sensitively  on  the  density  profiles  for  electrons  and  holes,  [Eq.  (19)] 
and  unless  these  densities  are  accurately  known,  it  is  Inconceivable  that  we 
shall  obtain  the  correct  sign  of  charge  on  EHD, 

From  Eq.  (21)  we  can  write  down  for  the  difference  (y^  - ye) ; 

y.  - u - y,  - y - 2A$  (22) 

n e n e 

T*'e  sign  and  magnitude  of  charge  is  completely  determined  by  the  difference 

y^  - ue«  Since  twice  the  dipole  layer  is  not  expected  to  remove  the  differ- 
ence y.  - y , the  aforementioned  conclusions  about  charge  will  remain  un- 
it e 

changed. 


II.  SELF-CONSISTENT  CALCULATION 

There  are  essentially  two  reasons  for  doing  a self-consistent  calculation 
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for  the  surface  properties  of  EHL.  In  view  of  the  reasonable  agreement  between 

the  variational  calculation  of  Vashlshta  et  al. and  the  experimental  meas- 

(39  41  42) 

urement  of  surface  tension  in  Ge(4;2)  * ’ , one  would  like  to  know  what  ef- 

fect the  higher  kinetic  and  exchange-correlation  gradient  corrections  have.  The 
other  motivation  for  doing  the  self-consistent  calculation  stems  from  the  fail- 
ure of  all  the  variational  calculations^®  ^ in  providing  a satisfactory  answer 
to  the  question  of  charge.  Let  us  remind  our  reader  that  these  variational  cal- 
culations yield  a positive  charge  on  EHD  in  Ge(4;2)  and  a negative  charge 

in  Ge(l»2)^^.  These  features  are  in  complete  contradiction  with  the  findings 

08) 

of  Pokrovsky  and  Svistunova'  ; they  measure  a negative  droplet  charge  in  Ge(4;2) 
and  a positive  charge  in  Ge(l;2).  Further,  on  the  basis  of  energetic  considera- 
tions we  showed  that  the  EHD  charge  would  be  negative  in  Ge(l;l),  whereas  the 
variational  calculation^^  yields  a positive  charge  on  the  EHD*  The  failure  of 
variational  calculations  is  due  to  the  fact  that  the  dipole  layer  is  overesti- 
mated. Considering  how  sensitive  dipole  layer  is  to  the  density  profiles  of 
electrons  and  holes,  the  way  to  improve  upon  the  variational  calculations  is  to 
obtain  these  densities  more  accurately.  Another  serious  shortcoming  in  varia- 
tional calculations  arises  because  of  truncation  of  gradient  expansion,  as  a re- 
sult of  which  one  misses  out  the  effect  of  Friedel  oscillations. 

All  these  shortcomings  can  be  remedied  by  the  elegant  self-consistent  scheme 
(44) 

of  Kohn  and  Sham  . It  allows  an  exact  treatment  of  noninteracting  kinetic 
energy,  and  consequently  one  does  not  miss  out  the  effect  of  Friedel  oscilla- 
tions arising  from  kinetic  energy.  The  limitation  of  the  self-consistent  pro- 
cedure lies  in  how  well  one  can  include  exchange  and  correlation  effects.  The 
Koin-Sham  formalism  is  based  on  the  following  three  equations: 

2 

(-  5 72  + Vjt[(n1>  ; r])  *J(r)  - *“(r)  (23) 
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d 


which  resembles  Che  Schrodlnger  equation.  The  first  term  on  the  l.h.s.  of 
Eq.  (23)  arises  from  the  noninteracting  kinetic  energy  of  component  of  EHL. 
The  second  term  on  the  l.h.s.  is  the  effective  potential  felt  by  a particle  of 
£th  component.  It  depends  on  the  densities  of  all  components,  {n^},  and  can  be 
written  as 


Vt({n1>:  r] 


+ 


«nt(r) 


(24) 


where  C^(r)  is  the  electrostatic  potential  of  l6*1  component,  given  by  Eq.  (12), 
and  6E  /6n^(r)  is  the  contribution  arising  from  the  exchange  and  correlation 
energy  of  EHL. 

The  third  equation  constitutes  an  expression  for  the  density  of  com- 
ponent in  terms  of  eigensolutlons, ^(r),  of  Eq.  (23),  and  it  reads 

nt(r)  - E |*“(r)|2  (25) 


where  a denotes  the  lowest  occupied  eigenstates.  Equations  (23)-(25)  consti- 
tute a set  of  self-consistent  equations  for  the  tc^  component.  Their  solu- 
tion yields  the  densities  (n^(r)},  which  are  the  basic  quantities  in  density 
functional  formalism.  It  is  apparent  from  Eqs.  (23)-(25)  that  these  are  a cou- 
pled set  of  equations  - the  potential  of  component  depends  not  only  on  the 
density  of  l**1  component  but  also  on  the  densities  of  all  other  components. 

This  is  the  primary  source  of  difficulty  in  solving  Eqs.  (23)-(25)  for  each 
component  of  EHL. 

It  is  evident  from  Eqs.  (23)— (25)  that  the  central  quantity  in  the  self- 
consistent  method  is  the  exchange-correlation  potential,  6Ejcc[{ni)]/6nl(£) . 
Assuming  that  the  densities  vary  slowly  in  space,  we  may  approximate  Ex(.  by 
the  local  density  term,  whereby 


E 


xc 


di{,.(t”i» 


(26) 
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The  quantity,  ({n^}),  is  the  total  exchange-correlation  energy  density.  It 

is  obtained  from  the  ground  state  energy  calculation  for  a system  of  uniform  den- 
sities {n^},  and  then  each  ii^  is  replaced  by  a varying  density  n^(r).  Eq.  (26) 
is  commonly  known  as  local  density  approximation  (LDA).  The  appropriateness 

of  LDA  for  the  determination  of  ground  state  energy  was  first  discussed  by  Kohn 
(44)  (53) 

and  Sham  . Tong  and  Sham  showed  that  in  LDA  the  density  n(r)  for  atoms 

agrees  well  with  the  Hartree-Fock  calculations.  Using  an  exact  sum  rule,  Vanni- 
(54) 

menus  and  Budd  have  argued  strongly  for  the  correctness  of  n(rj  obtained  by 

Lang  and  Kohn  in  LDA^^.  Thus,  it  seems  that  the  local  density  approximation 

is  indeed  a good  approximation  for  the  determination  of  density. 

Functional  differentiation  of  Eq.  (26)  leads  to  the  following  expression 

xc 

for  exchange  and  correlation  potential,  } ) : 


36  (( n . } ) 

({ni})  " 3^0— 


(27) 


Since  the  particles  belonging  to  component  may  reside  in  an  anisotropic 
elliptic  band,  the  kinetic  energy  term  in  Eq.  (23)  may  have  the  form 


2 2 2 2 

iL  = _L_  (11  + II  > + — i — i_ 

2m  - 2m,  „ \ 2 . V 2m,  . . 2 * 

i t,t  3x  3y  l, l 9z 


(28) 


where  m.  and  m are,  respectively,  the  tranverse  and  longitudinal  masses  of 
C *.  y A> 

a particle  belonging  to  the  component.  Under  the  volume  conserving  trans- 


formation 


where 


- 1 1 I 

x ■ y ^ x ; y - y ^ y ; z - y 3 z 


(29) 
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Equation  (28)  takes  up  the  form 


V_ 

2a. 


2S„ 


(30) 


where  the  kinetic  energy  operator,  i8  given  by 


2 2 

v2  = i-  +11  + 1_ 

- - 2 ^ 2 2 » 


2 ' 2 ’ 2 
3x  3y  3g 


(31) 


and  m^,  the  density  of  states  mass  for  a particle  of  1^  component,  is  de- 
fined by 


£ ■ (m  m )3* 

v 2,,t 


(32) 


Henceforth,  we  shall  drop  the  tilde  from  the  kinetic  energy  operator, 

and  omit  the  bar  over  the  density  of  states  mass,  m^. 


Assuming  that  the  EHL  surface  is  perpendicular  to  the  z-axis,  the  expres- 
sion for  ij»°(r)  deep  inside  the  liquid  is 


*“(r)  « sin(kz  - yOO)  e1(kxx  + kyy) 


(33) 


where  k^,  k^,  k are  the  three  components  of  wave  vector  !k  and  r(k)  is  the 
phase  shift.  Substituting  Eq.  (33)  in  Eq.  (23),  we  get 


€i  " \ (kx  + kJ  + fc2)  + V{5i}  ; ‘ “] 


Using  Eq.  (24)  in  Eq.  (34)  we  obtain 

.2 


**  9 2 2 f{n,  }J 

El  " 2^  (kx  + ky  + k > + V(-  m)  + —irt 


(34) 


(35) 


Substituting  Eq.  (27)  in  Eq.  (35),  we  find 


■ 3,  <k^  + ky  * k2>  + h*<-  •>  + r<«i» 


(36) 
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where  is  the  exchange-correlation  contribution  to  the  bulk  chemical  poten- 
dal  of  1th  component.  Eq.  (36)  holds  only  when  the  wave  function  is  given  by 
Eq.  (33).  In  general,  i|i°(r)  has  the  form 

o.  . . . i^x  + kyy)  • (37) 

- 4»t(k,z)  e ^ ™ 


Substituting  Eqs.  (36)  and  (37)  in  Eq.  (23)  and  expressing  length  in  excitonic 

v4 

Bohr  radius,  the  energy  in  excitonic  Rydberg  (1  Rydberg  “ j 2^  we 

2ft  K 

(-  ^ ^ + Viff[{ni}  ; 2l)  (k»  z)  - ^ (k2  - kjt>  ^(k,z)  (38) 

nl  dz  *■ 

• ^ 2 

where  - 3w  n0/v  , and  m is  the  reduced  mass.  For  the  effective  potential 
F*  * ’ r 

of  the  component,  we  obtain  the  expression. 


Vf[(ni>  ; z]  - 8*  C 


, r - 


(z  - z')  [n  (z' ) - n.  (z')] 


(39) 


- Ul({H1»  + wJC({nt»  , 


where  n#(z)  and  n^(z)  *re,  respectively,  the  electron  and  hole  densities,  and 
P is  the  bulk  chemical  potential  of  tth  component.  The  exchange-correlation 

JL 

potential,  u*c,  is  given  by  Eq.  (27). 

From  the  ground  state  energy  calculations  for  a homogeneous  EHL,  we  know 
that  to  a very  good  approxisutlon  the  electrons  and  holes  contribute  equally 
to  exchange-correlation  energy ( 26,37 * . Therefore,  it  is  reasonable  to  divide 
exchange  and  correlation  contribution  in  the  bulk  equally  among  electrons  and 
holes.  In  that  case,  Eq.  (39)  simplifies  to 


Vjff [{n±>  : *] 


+ Uj  (nt(z)) 


(z  - z')  [ne(z’)  - n^z')]  - ^(nj) 

(AO) 


-18- 


The  exchange-correlation  potential  is  given  by 


± <v>i» 

*> 


where  e*C(n.)  is  the  exchange-correlation  energy  per  particle.  Defining  a 


local  value  of  rg  by 


n£(z)  - 3/4*  rg  £ (z) 


the  expression  for  e^  (n^)  that  we  have  used,  reads 


xc.  . 
cl  (ni) 


a Z ...  b(i)-2  ( . 

r-T-  - c(i)  r , , r .(z)  < r 

r ,(z)  1-1  «.1  3*4  c 

3 1 1 


0. 5*B 


r , (z)  + C + r (z)  * rs,t 

8)'  O 8j* 


, r ,(z)  > r 


In  Eq.  (43)  a is  related  to  the  coefficient  of  exchange  energy,  the  c(i)  are 


the  coefficients  of  polynomial  fit  to  correlation  energy,  and  Bq  and  Cq  are 


the  coefficients  of  Wlgner  fit  to  correlation  energy.  The  coefficients  of 


polynomial  and  Wlgner  fits  to  correlation  energy  are  obtained  from  fully  self- 


consistent  calculation  of  Vashlshta  et  al.  . Their  calculation  includes 


the  effects  of  anisotropy  of  conduction  bands  as  well  as  multiple  scattering 


of  e - e,  e - h and  h - h.  The  binding  energy,  equilibrium  density,  the  cri- 


lcal  temperature,  obtained  by  Vashlshta  et  al.,  agree  very  well  with  experi- 


mental results.  Values  of  a,  c(l),  Bq  and  CQ  are  given  in  Tables  III  and  IV. 


Performing  the  sum  over  bands  and  the  lowest  states  in  Eq.  (25)  [see  Ap- 


pendix I]  ve  obtain. 


v .2  2 2 

t(Z)  - dk(kFt  - kZ)  [*t(k,Z)r  . 

IT  ' 0 


The  EHL  in  Ge(4;2),  Ge(l;2),  Sl(6;2)  and  Si(2;2)  consists  of  three  dis- 


F 


■ 


tinct  componenta  - electrons,  light  holes  and  heavy  holes.  Corresponding  to 
each  component  we  have  a set  of  equations  like  Eqs.  (38),  (AO)  and  (44).  It 
is  well  known  that  density  of  states  mass  of  light  holes  is  much  smaller  than 
that  of  heavy  holes  with  the  consequence  that  the  light  hole  contribution  to 
equilibrium  density  is  much  smaller  than  that  of  heavy  holes.  Therefore,  it 
is  reasonable  to  solve  only  the  set  of  equations  for  electrons  and  heavy  holes. 
In  the  set  of  equations  for  holes  we  include  the  kinetic  contribution  of  light 
holes  by  defining  an  effective  hole  mass,  m^: 


m u 3/2 

+ '2'3  • 


(43) 


where  m ^ and  m^  are,  respectively,  the  light  and  heavy  hole  masses  in  Ge 
and  Si.  The  use  of  m^  in  Eq.  (38),  for  example.  Insures  an  exact  treatment  of 
kinetic  energy  of  holes  in  a homogeneous  EHL.  In  Ge(l;l)  and  Sl(2;l),  there 
is  only  one  hole  band,  with  the  consequence  that  there  are  only  two  sets  of 
coupled  equations  to  be  solved. 

We  start  with  exponential  trial  density  profiles  for  electrons  and  holes, 
and  construct  their  effective  potentials  from  Eq.  (40).  Having  obtained  the 
effective  potentials,  we  solve  Eq.  (38)  for  i^(k,z)  of  each  component  and  then 
use  them  in  Eq.  (44)  to  construct  new  densities  (n^(z)}.  This  procedure  is 
continued  until  (nA(z)}  converge  satisfactorily.  The  degree  of  convergence 
that  we  attain  for  each  density  profile  is  better  than  0.5Z  of  the  mean  density. 
A further  check  on  convergence  is  provided  by  the  one-dimensional  analog  of 
Friedel  sum  rule^"^  [see  Appendix  I]  which  is  very  well  satisfied  in  our  cal- 
culations [for  details  of  our  numerical  procedure,  see  Appendix  II]. 

The  self-consistent  density  profiles  in  Ge(4;2),  Ge(l;2),  Ge(l;l),  Sl(6;2), 
Si ( 25  2) , and  Si(2;l)  are  shown  in  Figures  (1)  - (6).  Tables  V and  VI  contain 
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normalized  values  of  n (z ) and  n,  (z ) in  the  six  systems.  It  is  evident  from 

e n 

Figures  (1),  (3),  (4)  and  (6)  that  the  hole  density  spills  out  more  than  the 
electron  density  in  Ge(4;2),  Ge(l;l),  Si(6;2),  and  Si(2;l).  The  excess  spill- 
ing of  holes  is  due  to  the  fact  that  in  these  four  systems  the  binding  energy 
of  holes  is  less  than  that  of  electrons.  It  further  indicates  that  in  Ge(4;2), 
Ge(l;l),  Si(6;2)  and  Si (2; 1)  one  would  expect  the  EHD  to  carry  a negative 
charge.  In  Ge(l;2)  and  Si(2;2),  there  is  an  excess  leakage  of  electrons  than 
holes,  indicating  that  the  holes  have  a greater  binding  energy  than  electrons. 
Consequently,  the  EHD  is  expected  to  sustain  a positive  charge  in  Ge(l;2)  and 
Si(2;  2). 

It  is  apparent  from  Figs.  (l)-(6)  that  the  larger  the  binding  energy  of  a 
component,  the  bigger  the  amplitude  of  the  first  peak  in  Friedel  oscillations. 
This  trend  is  in  line  with  the  results  of  Lang  and  Kohn  for  a metal  surface^^ 
We  would  like  to  mention  here  that  the  Friedel  oscillations  in  effective  poten- 
tial are  not  as  pronounced  as  in  electrostatic  potential. 

III.  SURFACE  ENERGY 

The  surface  tension  of  EHL  consists  of  three  contributions ; kinetic  (o  ) , 

s 

electrostatic  (°eg)»  an<*  exchange  and  correlation  (o^). 

o-a  + o + o (46) 

s es  xc 

Each  of  these  terms  is  further  made  up  of  contributions  arising  from  elections 

and  holes.  Using  the  superscript  t to  denote  a component  of  EHL,  we  shall  give 

it  i 

below  the  expressions  for  o , o and  o 

s es  xc 

The  expression  for  a*  is  obtained  in  the  same  way  Lang  and  Kohn  did ^55  \ 

s 

In  terms  of  phase  shifts  and  effective  potential,  kinetic  contribution  of  i^ 
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component  to  surface  tension  is 

p ^pmr  rkFC  —2  2 7T 

°8[nji] " ~t~jl  dk  k(kFrk  ){f " Y(k)} 

11  md  ° 

00 

- dz(V®ff[{n1};z]-V®ff[{ni};-)]  )n£(z)  (47 

In  evaluating  it  is  extremely  important  that  the  Friedel  sum  rule^^  [see 
Appendix  I ] be  well  satisfied.  Violation  of  Friedel  sum  rule  can  lead  to 
serious  errors  in  the  kinetic  contribution  to  surface  energy. 

The  exchange  and  correlation  contribution  of  component  is  given  by 


o [nj 
xc  l 


|°°  dz  n 


a(z)[^C(ni)-£-(na)]  , 


where  (n^)  is  obtained  from  Eq.  (43) . 

The  electrostatic  contribution  to  surface  tension  can  be  written  as 


where 


Ogg  [(ni>]  “ -j  j”  dz  <p(z)  [ne(z)-nh(z)  ] 

—OO 

4>(z)  - 8 it  |*dz' (z-z')  [ne(z')-nh(z')  ] 


Using  the  self-consistent  density  profiles,  effective  potentials  and  phase 
shifts  we  evaluate  the  electrostatic, kinetic, and  exchange-correlation  contri- 
butions to  surface  tension.  Table  VII  contains  values  of  surface  tenstion,  a, 

for  EHL  in  Ge(4;2),  Ge(l;2),  Ge(l;l),  Si(6;2),  Si(2;2)  and  Si(2;l).  We  find 
-4  2 

that  o » 3.7  x 10  erg/cm  in  Ge(4;2).  It  agrees  favorably  with  the  measure- 

C39}  (41)  . (42) 

ments  of  Westervelt  et  al.  ' , Etienne  et  al.  and  Staehli  . Note,  the 
self-consistent  value  of  o in  Ge(4;2)  differs  appreciably  (almost  a factor  of 
4)  from  the  values  obtained  by  Rice^^  , Sander  et  al^^\  and  Reinecke  and 
Ying^®^.  It  substantiates  similar  claims  made  by  Lang  and  Kohn  in  the  context 


rr— — 

of  a metal  surface^^.  However,  when  one  includes  in  the  variational  calcu- 
lation the  gradient  correction  to  exchange  and  correlation  energy (50»51) ^ tfie 
value  of  o in  Ge(4;2)  is  found  to  be  3.5  x 10  erg/cm  . Such  an  agreement 
between  self-consistent  and  the  variational  calculations  is  rather  fortuitous. 
Within  the  framework  of  density  functional  formalism,  Rose  and  Shore  have  car- 
ried out  a partial  self-consistent  calculation  for  EHL  surface  in  Ge(4;2); 
the  authors  report  a value  of  2.6  x 10  ^ erg/cm^  for  EHL  surface  tension^^. 

We  find  that  the  value  of  o in  Ge(l;l)  is  a factor  of  twenty  smaller  than  in 
Ge(4;2). 

-4 

In  Si (6; 2),  the  surface  tension  of  EHL  is  calculated  to  be  87.4  x 10 

erg/cm  v . It  is  again  a factor  of  three  larger  than  the  variational  re- 

(31) 

suit  of  Reinecke  and  Ying  . The  self-consistent  result  for  a in  Si(2;l) 

is  a factor  of  eight  smaller  than  in  Si(6;2).  Considering  the  important  role 

played  by  surface  tension  in  the  determination  of  coexistence  curve,  kinetics 
of  EHD  formation,  and  in  establishing  the  validity  of  density  functional  form- 
alism, there  is  a dire  need  for  conclusive  experiments  on  EHL  surface  in  the 
various  configurations  of  Ge  and  Si. 

IV.  CHARGE  ON  ELECTRON-HOLE  DROP 

Calculations  for  the  ground  state  energy  of  a homogeneous  EHL  lead  us  to 
believe  that  the  constituents  of  EHL,  namely,  electrons  and  holes,  have  differ- 

j 

(20-26) 

ent  binding  energies  v ’ . In  the  presence  of  a surface  the  difference  in 

binding  energies  results  in  dissimilar  density  profiles  for  electrons  and  holes, 

and  consequently  in  a dipole  layer.  The  sign  of  dipole  layer  is  dictated  by 

the  energetics  of  electrons  and  holes,  since  the  purpose  of  dipole  layer  is  to 

(37) 

reduce  the  difference  in  their  binding  energies 
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Following  Che  conventional  definition  of  dipole  layer,  we  write 


A$  --8w 


j" dz  z[ne(z) 


nh(z))]  Ry  . 


It  is  implicit  in  Eq.  (50)  that  the  total  number  of  electrons  is  equal  to  the 
total  number  of  holes.  Using  the  self-consistent  electron  and  hole  densities, 
we  calculate  from  Eq.  (50)  the  values  of  dipole  layer  in  the  six  configurations 
of  Ge  and  Si.  These  values  are  given  in  Table  VII.  Evidently,  holes  have  less- 
er binding  energy  than  electrons  in  Ge(4;2)  and  Si(6;2).  In  order  that  the  dif- 
ference between  the  binding  energies  of  holes  and  electrons  be  reduced  by  the 
presence  of  EHL  surface  the  sign  of  dipole  layer  should  be  positive  in  Ge(4; 2) 
and  Si(6; 2) , which  is  indeed  the  case.  Similar  considerations  lead  us  to  con- 
clude that  A$  should  be  negative  in  Ge(l; 2)  and  Si(2; 2)  and  positive  in  Ge(l; 1) 
and  Si(2; 1) . 

Charge  on  EHD  is  determined  by  the  difference  in  the  chemical  potentials 
of  electrons  and  holes.  In  terms  of  bulk  chemical  potentials  and  dipole  layer, 
this  difference  is  given  by  (see  Eq.  (22)), 


y,  - y - y - y - 2A$  (51) 

h e n e 

It  is  apparent  from  Table  VII  that  the  difference  y.  - y is  positive  in  Ge(4:2) 

n e 

and  Si(6;2),  which  means  that  the  holes  are  less  tightly  bound  than  electrons. 

At  a finite  temperature  (T  < T , the  transition  temperature  for  EHL)  there  will 

be  an  excess  thermionic  emission  of  holes,  and  as  a result  the  electron-hole 

droplet  will  acquire  a negative  charge  whose  magnitude  will  be  determined  by 

the  condition  that  at  equilibrium  the  work  functions  for  electrons  and  holes  be- 
(52  58) 

come  equal  * . A similar  situation  occurs  in  Ge(l;  1)  and  Si(2;  1),  wherein 

the  EHD  develops  a negative  charge.  In  Ge(l;2),  however,  we  find  that  the  dlf- 

(52) 

ference  u.  - y is  negative,  implying  that  the  EHD  sustains  a positive  charge  , 
n e 
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whereas  in  Sl(2;2)  we  find  that  Che  difference  u,  - v is  extremely  small. 

h e 

Within  the  limits  of  accuracy  of  our  calculations,  we  obtain  a neutral  EKD  in 

Si(2  52) . Insofar  as  the  sign  of  charge  is  concerned,  our  results  in  Ge(4;2) 

(38) 

and  Ge(l;2)  agree  with  the  experiment  of  Pokrovsky  and  Svi3tunovav  . Recent 

(591 

measurements  of  Nakamura  has  also  shown  that  EHD  is  negative  in  Ge(4;2)v  No 

measurement  of  charge  has  yet  been  reported  in  Ge(l;l)  or  any  of  the  configu- 
rations in  Si. 

« 

There  are  two  idealized  situations  in  which  one  can  determine  the  magnitude 
of  charge  on  EHD.  First,  under  pulsed  excitation,  i.e.,  when  the  medium  outside 
of  EHD  is  regarded  as  vacuum,  the  magnitude  of  charge,  Q,  is  calculated  from  the 
equation 


Q 


[/2  + R|uh  - ufi|  - ll2 


(52) 


where  all  the  quantities  are  measured  in  reduced  units.  For  a typical  drop  of 

-4  i i 

radius  R - 284  a^(-  5 x 10  cm),  we  find  that  the  EHD  carries  a charge  of  - 2 7 1 e j 

in  Ge(4!2)^^.  Variational  calculations  of  Rice^2\  and  Reinecke  and  Ying^8^ 

yield  - 18|e|  and  + 6|e|,  respectively.  Calculation  of  Rose  and  Shore  also 

yields  a negative  charge  on  EHD  in  Ge(4;2)^2\  The  authors,  however,  do  not 

quote  the  magnitude  of  charge. 

The  second  situation  occurs  when  the  EHD  becomes  surrounded  by  ionized 
carriers.  The  problem  is  now  much  more  complicated  because  of  screening  by 
free  carriers.  Assuming  that  the  screening  length  is  much  less  than  the  radi- 
us of  EHD  and  treating  ionized  carriers  as  a classical  gas.  Rice  obtains  the 

(37) 

following  expression  for  the  droplet  charge  : 
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where  Debye-Hiickel  screening  length  is  given  by 

>5 

.V  . 

r 


kj»T 

r-JL  > 
',8trni/ 


(54) 


All  the  quantities  in  Eqs.  (53)  and  (54)  are  measured  in  reduced  units.  Under 

conditions  of  thermodynamic  equilibrium  at  4.2*K,  the  density  of  free  carriers, 
•12  “3 

n^  ■ 2 x 10  cm  . Using  this  value  of  n^  we  obtain  a droplet  charge  of  - 700 | e | 
(52) 

in  Ge(4;2)  . Measurement  of  Pokrovsky  and  Svistunova 

yields  a value  of  - 100 1 e | ^ 32  ^ , while  Nakamura  obtains 

(S9) 


a value  of  -400|e|  at  1.6°K. 


The  apparent  discrepancy  be- 


tween theory  and  experiment  is  mainly  due  to  the  fact  that  theoretical  estimates 
of  charge (-^,52)  are  M<je  un<jer  the  simplifying  assumption  that  the  screening 
length  is  much  less  than  the  radius  of  the  drop.  Also,  the  real  experimental  con- 
dition may  not  conform  well  to  the  above  mentioned  assumption.  It  should  also 

be  borne  in  mind  that  in  the  experimental  determination  of  charge,  there  are  am- 

(59) 

blguities  arising  from  the  assumption  of  collision  time 

It  is  evident  from  Eqs.  (52)  and  (53)  that  in  order  to  obtain  the  magnitude 
of  charge  one  needs  to  know  the  radius  of  EHD.  Experimental  measurements  of 
droplet  radius  have  not  yet  been  made  in  Ge(l;2),  Ge(l;l),  Si(2;2)  and  Si(2;l), 
which  prevents  us  from  estimating  the  charge  on  EHD  in  these  systems. 
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V 


CONCLUSION 


We  have  shown  that  in  order  to  resolve  the  question  of  EHD  charge  it  is  im- 
portant to  treat  the  kinetic  energy  of  electrons  and  holes  properly.  Since  the 
variational  method  provides  a poor  approximation  for  kinetic  energy ^87’48, 80 ^ , 

it  fails  to  account  for  the  experimental  observation  that  the  EHD  sustains  a 

(38) 

negative  charge  in  Ge(4;2)  and  a positive  charge  in  Ge(l;2)v  . By  solving  a 
set  of  coupled  Kohn-Sham  equations  for  the  constituents  of  EHD,  we  indeed  find 
that  the  EHD  is  negatively  charged  in  Ge(4;2)  and  positive  in  Ge(l;2).  Although 
the  Kohn-Sham  procedure  Involves  much  more  numerical  work  than  the  variational 
method,  the  former  has  the  advantage  that  it  treats  kinetic  energy  exactly  and 
also  includes  the  effect  of  Friedel  oscillations.  On  the  basis  of  self-consist- 
ent calculation  we  predict  that  the  EHD  carries  a negative  charge  in  Ge(l;l). 

Similar  features  are  also  expected  to  occur  in  silicon  under  a uniform  stress 
along  <100>  direction.  In  Sl(6;2)  and  Si(2;l)  our  calculations  show  that  the  EHD 
sustains  a negative  charge.  However,  in  Si(2;2)  we  find  that  the  electron-hole 
drop  is  almost  neutral. 

As  regards  surface  tension  of  EHD,  the  value  in  Ge(4;2)  agrees  favorably  with 

the  experimental  estimates  of  Westervelt  et  al.  Etienne  et  and  of 

Staehli^tae  self-consistent  estimate  of  a in  Ge(l;l)  is  found  to  be  a factor 

of  twenty  smaller  than  in  Ge(4;2).  The  surface  tension  in  Si(6;2)  is  calculated 
-4  2 

to  be  87.4  x 10  erg/cm  , a factor  of  eight  larger  than  in  Si(2;l).  The  values 

-4  -4,2 

of  o in  Ge(l;2)  and  Sl(2;2)  are  found  to  be  1.0  x 10  and  32.8  x 10  erg/cm  , re- 
spectively. 

One  of  the  approximations  in  this  work  amounts  to  dividing  equally  between 
electrons  and  holes  the  exchange  and  correlation  energy  per  e - h pair  in  the 
bulk  of  electron-hole  liquid.  This  is  a reasonable  approximation,  since  we 
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know  from  Che  ground  state  energy  calculations  for  a uniform  EHL  that  the  elec- 
trons and  holes  contribute  almost  equally  to  the  exchange-correlation  energy^^'-^) 
One  would  expect  it  to  be  the  case,  because  exchange  and  correlation  effects  a- 
rise  from  Pauli  principle  and  Coulomb  interaction.  To  avoid  confusion,  we  would 
like  to  emphasize  that  this  approximation  does  not  imply  that  the  exchange  and 
correlation  potentials  for  electrons  and  holes  are  equal  in  the  surface  region  of 
EHL. 

The  major  approximation  in  our  calculation  is  that  we  have  retained  only  the 
local  density  contribution  to  exchange-correlation  potentials  of  electrons  and 
holes.  It  may  seem  feasible  to  improve  the  local  density  approximation  by  includ- 
ing the  exchange-correlation  gradient  corrections.  The  first  gradient  correction 
to  exchange-correlation  energy  is  typically  of  the  form  j dr_  | Vn (_r ) ( /n^^(jr). 

It  can  be  easily  seen  that  the  potential,  obtained  by  taking  the  derivative  with 
respect  to  n(r),  behaves  pathologically  in  the  tail  region  of  density.  For  this 
reason,  it  is  not  possible  to  include  in  the  self-consistent  calculation  the  first 
gradient  correction  to  exchange-correlation  potential.  One  should  note,  however, 


that  such  a difficulty  does  not  occur  in  variational  calculations  that  include 
the  first  gradient  correction  to  exchange-correlation  energy,  which  is  well  be- 
haved everywhere.  Kohn  and  Sham  have  suggested  including  in  the  LDA  not  only 

the  first  gradient  correction,  but  also  the  higher  terms  in  the  gradient  expan- 
(44) 

sion  . Such  a calculation  will  involve  enormous  numerical  work.  Many  calcu- 
lations have  demonstrated  the  appropriateness  of  local  density  approximation  (LDA) . 
Gunnarson  et  al.  have  used  LDA  in  calculating  the  binding  energy  of  hydrogen 
atom^^.  Their  result  for  the  binding  energy  differs  from  the  exact  value  by 
less  than  two  per  cent.  Tong  and  Sham  have  shown  that  the  density  of 

atoms,  obtained  in  the  local  density  approximation,  agrees  well  with  the  Hartree- 
(53) 

Fock  calculations  . The  classic  work  of  Lang  and  Kohn  also  shows  that  the 
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LDA  is  a good  approximation  in  the  study  of  metal  surfaces 


(55) 


Recently,  Ando 


has  used  It  in  the  calculation  of  intersubband  separation  in  Si  Inversion  Layer 
His  results  are  in  excellent  agreement  with  experiment. 

In  view  of  the  fact  that  the  EHD  surface  is  a genuine  inhomogeneous  plasma  of 
electrons  and  holes,  it  is  an  ideal  system  for  testing  the  validity  of  density 
functional  formalism;  in  particular,  of  local  density  approximation  in  the  Kohn- 
Sham  procedure.  Thus  there  is  a dire  need  of  doing  careful  measurements  of  sur- 
face tension  and  charge  on  EHD  in  the  aforementioned  systems. 


(61) 
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APPENDIX  I 


Let  us  assume  that  the  EHL  has  dimensions  (L  , L , 2L  ) along  the  x,y  and  z 

x y z 

directions.  Further,  ve  suppose  chat  L , L >>  2L  and  that  Che  geometrical  sur- 

x y z 

face  of  EHL  lies  at  z ■ 0. 

The  Schrodinger-like  equation  for  the  1th  component  is  given  by 


*2  a2 


2 2 
1Tk 


{-  s.  7"2  + w ; z]}  vk’z)  * im~  vk*z) 


(A-l) 


1 3z 


For  Che  present  discussion  it  is  not  necessary  to  specify  the  form  of  potential 


V 

Multiplying  Eq.  (A-l)  with  ^£(k!,z),  interchanging  prime  and  unprimed  quan- 
tities in  the  last  equation  and  Chen  subtracting  the  resulting  equation  we  ob- 
tain. 


a 3*f(k,z)  3*,(k\z) 

az  [Vk’’z)  — a! Vk-z)  — a7 1 

- (k'2-  k2)i|»j(k'  ,z)4>^(k,z) 


(A-2) 


Integrating  Eq.  (A-2)  from  z to  • and  taking  the  limit  k'-*  k,  we  find 


f 21 

j a,-  iv.'’)] 


a*,(k,z)  ag»  (k,z)  a 4>,(k,z) 

-Si—  - bir  » <*-3> 


Deep  inside  the  liquid,  the  wave  function  takes  the  form. 


where 


^(k.z)  - At(k)  sin(kz  - y(k)) 


. mir  y(k)  . . 

k » ^ , m - even  Integer 

z z 


(A-4) 


and  0 < y(k)  < ir/2 


Taking  the  limit  z ■ - L^  in  Eq.  (A-3)  and  making  use  of  Eq.  (A-4)  we  get 


r 2 i 

j dz'  (^(k.z')]  I-  k(Lz  + y'(k))  ♦ j sin  2(kLz  + y 

* L 


(k))} 

(A-5) 
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where  y'(k)  = dy(k)/dk. 

Making  use  of  che  fact  that  the  wave  function  Is  either  symmetric  or  anti- 
symmetric about  z ■ - Lz>  l.e., 

-Lz 

dz  (k,  z)  ] 2 - | dz  (^(k.z)]2  - re— £ — , (A-6) 

L Llz  x y 

we  obtain  from  Eq.  (A-5)  the  expression  for  the  normalization  constant  A (k) : 


{ *t(k)}2  + .2  , 


(A- 7) 


where  ■ 2L^L^L^  is  the  volume  of  the  electron-hole  liquid. 

The  expression  for  the  density  of  component  [Eq.  (25)1  reads 

n,(z)  - 2v,  E I E | I2  . 
x y 


(A-8) 


where  the  factor  of  two  arises  from  spin  summation  and  comes  from  the  sum  over 

the  bands  acccmmodatire  the  particles  of  component.  The  lowest  filled  states 

are  characterized  by  m,  m and  m , such  that 
J x y 


(k  , k , k) 
x y 


2irm  2xm  . . 

( x L m7t  _ UJLL) 

V T » T • 2L  t / 


(A- 9) 


x y z z 

Now,  the  states  In  wave  number  space  lie  on  sheets  normal  to  k-  axis  and  are 

spread  according  to  Eq.  (A-9).  Since  L , L » 2L  , the  states  on  any  sheet 

x y z 

are  much  more  densely  packed  than  the  spacing  between  the  sheets.  It  Is  clear 

2 

from  Eq.  (A-9)  that  on  a given  sheet  there  are  (L  L /4ir  ) states  per  unit  area. 

x y 

We  take  the  occupied  states  to  lie  within  a hemisphere  of  radius  k^  in  the 


k > 0 half  space.  The  hemisphere  cuts  the  k-  axis  between  Rth  and  (R  + l) 
sheets  such  that 


th 


Si  - Sip  " Si  + 1 


(A- 10) 
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The  radius  of  che  circular  sheet,  which  accommodates  particles  of  mc^  state,  Is 
given  by 

•>  m.  o 1, 

(A-ll) 


<k?P  - <1  - W ■ 

z z 


Therefore,  the  expression  for  n^f^)  becomes, 


L L 


nt(z)  " 1 |At(k)|"l*t(k,Z)|2aJ  * kZ)  (A-12) 


m " * 

Changing  the  sum  over  m In  Eq.  (A-12)  Into  an  Integration  over  k and  using 
Eq.  (A- 7)  we  obtain 

v.  rk. 


nt(z)  - | ip  dk  (k2p  - k“) ^(k.z) j2  (A-13) 


where 


ib  (k,z)  -*■  T sm(kz  - y(k ) ) 

* £ z -*■  - L 

z 


Now,  total  ntimber  of  particles  of  component  Is  given  by 


(A-14) 


N-2v.EE  E 1 
l Imam 

x y 

Carrying  out  the  summations  In  Eq.  (A-15)  we  find 


(A-15) 


v k,  , 

N,  - 0.  (“M*)  1 1 “ T~ 
1 d 3n2  k3  L 


IP  z 


However , 


\ . Vn 

nd  3ir2 


Comparing  Eqs.  (A-16)  and  (A-17),  we  obtain 

*Ft 

"Fl“z  0 


Klp 


dkkf^--Y(k)}]  (A-16) 


(A-17) 


kIp  ■ Ef.u  + I"1'  dk  k(?  - '■<k)H  + 0(7-> 


(A-18) 


In  the  limit  z ■*  - L^,  we  obtain  from  Eq.  (A-13) 
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(A-19) 


r v 


v.k 


na(-Lz)  ” “M2  + OC-V-) 


3ir 


Owing  to  che  face  chac  no  free  charges  can  exisc  deep  inside  Che  EHL,  ve  de- 


mand 


N, 

--L  ‘57 

z d 


and  Cherefore  from  che  lasc  Cwo  equations  we  find 


Vfi 

3.2 


(A-20) 


% ’ St  + <A-21> 

In  order  chac  Eq.  (A-21)  be  conslscenC  wich  Eq.  (A-18) , we  musC  demand  thac 
for  each  componenC  t, 

dk  kY(k)  - j kj£  (A-22) 

o 

and  chus 

[K,  _2  2 i 2 

nt(z)  - \ri  dk  (kjt  - kz)  [^(k.or  (A-23) 

* o 

Equation  (A-22)  is  the  one-dimensional  analog  of  Frledel  sum  rule^56^ . It 
provides  a check  on  the  phase  shift  obtained  from  che  wave  functions.  Further, 
by  satisfying  this  sum  rule  one  Insures  that  both  the  electron  and  hole  densi- 
ties are  equal  to  the  mean  density  deep  inside  the  liquid. 
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APPENDIX  II 


We  shall  deal  here  with  Che  numerical  procedure  used  in  solving  a sec  of  cou- 
pled Kohn-Sham  equations  for  Che  case  of  EHL  surface.  Our  procedure  Is  a gen- 
eralization of  that  developed  by  Lang  and  Kohn  for  a metal  surface^"*^. 

We  stare  with  trial  density  profiles  for  electrons  and  holes.  These  pro- 
files decay  exponentially  outside  the  surface  and  approach  the  mean  bulk  density 
well  inside  the  surface.  Further,  the  parameters  in  the  initial  profiles  are  so 
chosen  that  we  obtain  overall  charge  neutrality: 


r 


dz  n (z) 
e 


(B-l) 


Using  the  trial  profiles  in  Eq.  (40)  we  calculate  the  effective  potentials 
for  electrons  and  holes.  Substituting  the  effective  potential  in  Eq.  (38),  we 
solve  it  for  (k,z)  in  the  range  0 to  ic  . Imposing  the  condition  that  the 

*>  f fc 

vavefunctlon  should  approach  a sine  wave  deep  inside  the  electron-hole  liquid, 


1 . e . , 


we  obtain  y(k)  from  Eq. 
rule(56)  [Eq.  (A-22)], 


♦t(k»*)  g _ Afc(k)  sin(kz  - y(k))  (B-2) 

(B-2)  and  check  if  the  phase  shifts  satisfy  the  phase 

f"  * -f  «-3> 


0 

for  every  component.  If  the  phase  rule  is  not  satisfied  to  a few  percent,  we 
discard  the  initial  trial  profiles,  choose  a new  set  of  parameters  in  the  ini- 
tial profiles  and  follow  the  aforementioned  procedure  until  Eq.  (B-3)  is  almost 
satisfied.  Note,  imposition  of  phase  rule  is  another  way  of  Insuring  that  the 
densities  approach  the  correct  bulk  values  deep  Inside  the  EHL.  The  asymptotic 
form  for  wavef unctions  [Eq.  (B-2)]  is  ideally  achieved  only  if  z -*■  - •».  How- 
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ever,  in  practice  it  is  not  feasible  to  go  to  an  arbitrarily  long  distance  inside 

the  liquid.  Therefore,  one  has  to  Introduce  correction  terms  to  Eq.  (B-2).  In 
order  to  deal  effectively  with  the  asymptotic  region  we  observe  that  the  den- 
sity of  component  assumes  the  following  asymptotic  form: 


n*(z)z 


,n£n  + 


3 cos  2CkFJ£z  - y(kF£)) 


(2kFtZ)i 


1 + o(-y-) 


(B-4) 


With  this  form  for  the  densities,  the  asymptotic 
ponent  becomes 


Viff[{V  5 Z]z 


r -2 

+ Z Pi,i' 
1 l ' 


cos 


value  of  potential  of  i com- 
v(XFl,))/z2  + 0(l/z3) 


(B-5) 

The  coefficients  V , are  determined  by  fitting  the  effective  potentials 

L %L 

2 

to  Eq.  (B-5).  Using  this  asymptotic  form  for  effective  potential  we  find  1/z 

correction  term  to  the  asymptotic  value  of  wavefunctlon  ifi£(k,z),  namely, 

Eq.  (B-2).  Now  the  wavefunctlon  is  fitted  to  a form  given  in  Eq.  (B-2)  plus 
2 

the  1/z  correction  term.  From  such  a fit  we  extract  the  value  of  normaliza- 
tion constant  A£(k). 

The  expression  for  density  involves  i^(k,z),  which  must  attain  a sinusoid- 
al form  deep  inside  the  liquid.  In  order  that  the  density  of  a given  component 
attains  its  correct  bulk  value,  the  amplitude  of  the  sine  wave  must  be  unity. 
Thus,  the  knowledge  of  A£(k)  enables  one  to  construct  i|>£(k,z)  such  that 


i|>,(k,z)  - sin(kz  - y(k))  (B-6) 

We  substitute  these  4»£(k,z)  in  Eq.  (44)  to  obtain  the  density  n£(z).  When 

the  new  {n£(z)>  are  close  to  the  initial  density  profiles,  we  resort  to  a lin- 

0 

ear  response  procedure  to  bring  about  self-consistency.  Let  us  denote  by  ln£) 
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and  {n^  ' } the  trial  densities  and  the  densities  obtained  after  first  iteration. 


Then,  Che  aeC  of  Eqs.  (38),  (40),  and  (44)  [with  {n^}  replaced  by  {n^}  In  Eqs. 


(38)  and  (40)  and  (nj  j In  Eq. (44)]  may  be  caken  Co  define  Che  following  func- 


tlonal  relationship: 


{n^}  - F{(np  ; z] 


(B— 7) 


Using  Che  sec  of  denslCles  {n^  ' } , we  obcaln  from  Eqs.  (38),  (40),  and  (44) 


anocher  sec  of  denslCles  (n^  ' 


In  Che  splrlc  of  Eq.  (B-7)  we  wriCe 


{n4(2)}  - F[{n^1)}  ; z] 


(B-8) 


In  order  ChaC  {n^}  be  the  true  solutions  we  must  have 


(B-9) 


Let  us  assume  that  the  addition  of  {6n^}  to  {n^}  brings  about  self-con- 


sistency. Then  we  obtain  from  Eq.  (B-9): 


n<X)  + ; z] 


(B-10) 


Since  {n£X^}  are  close  to  the  true  solutions,  the  quantities{6n^X^}  must  be 


small  compared  to  {n^X^}.  Making  a Taylor  series  expansion  of  the  r.h.s.  of 


Eq.  (B-10),  we  find 


6F,[{n,(1)}  ; z] 


(1)  . ,U(D  _ (2)  , M f\  (1)  , 6FtItnl  } 

nt  + 6ni  n!  + i£i  j 6n!  (z  > 6n(i)(j?) 


dzf  (B-ll) 


We  choose  for  <Sn£  ' a linear  combination  of  derivatives  of  harmonic  oscill- 


ator functions: 


6ni[1)  ’ - JS  al(1)  Ql(z)  * 


(B-12) 


where  Q^(z)  are  derivatives  of  harmonic  oscillator  functions,  NH  the  number  of 


such  functions,  and  *^(1)  are  the  coefficients  whose  determination  shall 
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* 


enable  us  to  obtain  the  self-consistent  density  profiles.  In  order  that  (Q^(z)} 
form  a basis  set,  NH  must  be  infinite.  However,  in  practice,  one  needs  a reason- 
ably large  but  finite  number  (~20)  of  (^(z)  to  achieve  convergence.  Substituting 
Eq.  (B-12)  in  (B-ll)  and  using  the  definition  of  a functional  derivative,  we  ob- 
tain 


M NH  a ( j) 


n*  = i^1  "V-  (F[ni*ni-XQj’”1'F[nl*ni’"]+X,5UQj)  (B_13) 

where  X is  a small  parameter.  Multiplying  Eq.  (B-13)  with  harmonic  oscillator 
functions,  P^,(z),  and  integrating  over  z we  obtain 

r n\  M ^ a-»0)  r 

j dz  [n^1 ' (z)-n^  )(z)]?y(z)  m l l — j — J dz  , (z)  (F^  [n^.n^^-XQ^ , — ] 


- F[n1,njL,  — ]+X6i2Qj) 


(B-14) 


To  make  this  procedure  effective  one  has  to  choose  carefully  the  position 
and  width  of  harmonic  oscillator  functions.  It  is  our  experience  that  these 
functions  should  always  be  centered  near  the  surface  and  their  width  should  be 
on  the  order  of  surface  thickness. 

From  equations  such  as  (B-14)  we  obtain  the  coefficients  a^(j)  and  subse- 
quently use  them  in  Eq.  (B-12)  to  obtain  6n£^ . Straightforward  addition  of 
6n£X^  to  n£^  yields  the  self-consistent  values  of  density  of  component. 
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Table  1 - List  of  band  masses,  dielectric  constant  and  the  value  of  excitonic 
rydberg  in  Ge  (4;2)  and  Si  (6; 2),  and  ■ are,  respectively,  the  longi- 

tudinal and  transverse  masses  for  an  electron,  m . and  m,  the  optical  and 

oe  de  r 

density  of  states  masses  for  an  electron,  and  m^  and  m^  are  the  light 
and  heavy  hole  masses.  These  masses  are  given  in  units  of  bare  electron  mass. 
k is  the  dielectric  constant  of  the  system  and  Ex  is  the  value  of  excitonic 
tydberg.  The  values  of  band  masses,  k and  Ex  in  Ge  (1;2)  and  Si  (2;2)  are 
the  same  as  in  Ge  (4;2)  and  Si  (6;2),  respectively. 


SYSTEM 

"el 

■ * 

et 

■ 

oe 

a. 

de 

“lh 

"hh 

tc 

E (meV) 
X 

Ge  (4 ;2) 

1.58 

0.082 

0.120 

0.22 

0.042 

0.347 

15.36 

2.65 

Si  (6;2) 

0.9163 

0.1905 

0.2588 

0.32 

0.154 

0.523 

11.4 

12.85 

Table  II  - Values  of  constants  in  Ge  (1;1)  and  Si  (2;1).  The  masses  are 


measured  in  units  of  bare  electron  mass.  m^g  and  m^  are  the  density  of 

states  masses  for  an  electron  and  hole;  m and  m . are  their  optical 

oe  on 

masses,  m^  and  n>ht  are  the  longitudinal  and  transverse  masses  for  a hole. 
Ex  is  the  value  of  excitonic  rydberg  in  Ge  (1;1)  and  Si  (2;1). 


SYSTEM 

m , 
el 

m 

et 

m , 
de 

3 

O 

A 

mhl 

mht 

mdh 

m . 
oh 

Ex(meV) 

Ge(ljl) 

1.580 

0.082 

0.2198 

0.120 

0.040 

0.130 

0.088 

0.075 

2.65 

Si(2;l) 

0.9163 

0.1905 

0.3216 

0.2588 

0.1988 

0.2561 

0.2354 

0.2336 

12.85 
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Table  III  - Coefficients  of  exchange  and  correlation  energy  per  particle  in  unstressed 
Ge  and  Ge  under  uniform,  uaiaxial  stress  along  <111>  direction.  The  coeffi- 
cient of  exchange  energy,  a,  is  measured  in  excitonic  rydberg.  For  rs(Z)  < rQ 

correlation  energy  is  fitted  to  a polynomial  in  r (Z),  where  r (Z)  = {^-n(Z)} 

s so 

Coefficients  of  the  polynomial  fit  are  given  by  c(i)  = 0.5  * A(i)/b(i)  where 

b(.i)  = i + 0.25.  Correlation  energy  is  taken  to  be  of  the  Wigner  form 

for  r (Z)  > r . B and  C are  the  coefficients  of  the  Wigner  fit.  A(i)  and  B 
s o o o o 

are  expressed  in  excitonic  rydberg;  Cq  and  are  dimensionless. 


COEFFICIENTS 

Ge(4; 2) 

Ge(l;2) 

Ge(l;l) 

a 

-0.5681 

-0.7090 

-0.8297 

A(l) 

0.7212729 

0.4125300 

0.4683857 

A(2) 

4.1265679 

4.0192903 

1.8475582 

A(3) 

-9.6488301 

-9.4066370 

-1.5382064 

A(4) 

13.9998052 

13.2843140 

0.9232383 

A(  5) 

-10.3165219 

-9.2912027 

-0.2690568 

A(6) 

3.5668684 

3.0446329 

0.0282182 

A(7  ) 

-0.4567422 

-0.3725362 

0.0 

B 

-3.10247 

-3.27550 

-4.21484 

o 

C 

1.45138 

2.18922 

3.91350 

o 

r 

o 

2 

2 

3 

Table  IV  - Values  of  coefficients  of  exchange  and  correlation  energy  per  par- 
ticle in  Si(6;2),  Si(2;2),  and  Si(2;l).  The  coefficients  of  polynomial  fit  to 
correlation  energy  are  given  by  c(i)  = 0. 5*A(i)/b(i) , where  b(i)  = i + 0.25. 


Values  of  a. 

A(i) , B , C , and  r 
o o o 

are  tabulated  below. 

Note,  o,  A(i ) , and  B 

are  measured 

in  units  of  excitonic  rydberg.  C and  r 

o o 

are  dimensionless. 

COEFFICIENTS 

Si(6;2) 

Si(2;2) 

Si(2;l) 

a 

-0.5832 

-0.6891 

-0.8026 

A(l) 

0.7105611 

0.4000295 

0.4362952 

A(2) 

3.5159968 

4.8634251 

2.3587707 

A(3) 

-6.5141196 

-13.0590236 

-3.1982320 

A(4) 

8.0945300 

20.4316281 

3.4135213 

A(  5) 

-5.4580137 

-15.6070080 

-1.9050812 

A(6) 

1.7901528 

5.5015154 

0.4931726 

A(7) 

-0.2217641 

-0.7125849 

-0.0469333 

B 

o 

-3.32745 

-3.16211 

-4.21917 

C 

o 

1.76294 

1.95952 

3.93928 

r 

o 

2 

2 

3 

Table  V.  Self-consistent  density  profiles  of  electrons  and  holes  in 


Ge(4;2),  Ge(l;2),  and  Ge(l;l).  The  origin,  z = 0,  is  taken  to 

be  the  position  of  geometrical  surface  of  EHL.  The  distance, 

z,  is  measured  in  units  of  excitonic  Bohr  radius,  ax  (*177X)  , 

and  the  densities  in  units  of  a 3.  In  this  table  the  densities  have 
• x 

been  normalized  to  the  mean  e-h  pair  density.  Although  we  calculate 
the  density  profiles  at  points  separated  by  a distance  of  0.15  a^,  we 
quote  the  values  of  ng( z)  and  n^(z)  only  at  every  other  point. 

Ge(4-,2)  Ge(l;2)  Ge(l;l) 


z nfi(z)  n^(z)  ng(z)  n^(z)  ne(z)  n^(z) 


Table  V (Cont'd) 


Ge(4;2)  Ge(l;2)  Ge(l;l) 


z 

n (z) 
e 

nh<2) 

n (z) 
e 

nh(z) 

n (z) 

6 

nh(z) 

-0.75 

1.010 

0.992 

0.979 

0.999 

0.852 

0.797 

-0.45 

0.871 

0.859 

0.840 

0.864 

0.682 

0.649 

-0.15 

0.602 

0.610 

0.594 

0.606 

0.507 

0.500 

0.15 

0.315 

0.336 

0.344 

0.336 

0.344 

0.356 

0.45 

0.135 

0.139 

0.166 

0.150 

0.208 

0.248 

0.75 

0.046 

0.049 

0.069 

0.053 

0.122 

0.166 

1.05 

0.014 

0.017 

0.026 

0.017 

0.066 

0.106 

1.35 

0.002 

0.003 

0.008 

0.004 

0.027 

0.059 

1.65 

0.000 

0.001 

0.002 

0.001 

0.011 

0.031 

1.95 

0.000 

0.000 

0.001 

0.000 

0.006 

0.019 

2.25 

0.000 

0.000 

0.000 

0.000 

0.003 

0.014 

2.55 

0.000 

0.000 

0.000 

0.000 

0.001 

0.010 

2.85 

0.000 

0.000 

0.000 

0.000 

0.000 

0.007 

Table  VI  Self-consistent  electron  and  hole  densities  in  Si(6;2),  Si ( 2 ; 2 ) and 

O 

Si(2;l).  z is  measured  in  units  of  excitonic  Bohr  radius,  ax(=50A), 

-3  . 

and  the  densities  in  ax  . The  origin  is  taken  to  be  the  position  of 

the  geometrical  surface  of  EHL.  In  this  table  the  densities  have 

been  normalized  to  the  mean  e - h pair  density. 
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z 

Si(6;2) 

Si(2;2) 

Si(2 

;D 

n#(z) 

”h(l> 

ne(z) 

nh(2) 

n (s) 
e 

nh(z) 

-5.85 

0.996 

0.995 

1.006 

1.002 

1.011 

1.004 

-5.55 

0.992 

0.994 

1.000 

0.996 

1.017 

1.003 

-5.25 

0.987 

0.987 

0.994 

0.994 

1.014 

1.001 

-4.95 

0.994 

0.997 

0.992 

0.997 

1.010 

1.000 

-4.65 

1.004 

1.007 

0.992 

1.000 

1.002 

0.999 

-4.35 

1.012 

1.011 

1.000 

1.004 

0.990 

0.998 

-4.05 

1.013 

1.011 

1.008 

1.004 

0.979 

0.996 

-3.75 

1.016 

1.011 

1.013 

1.003 

0.974 

0.993 

-3.45 

1.002 

1.002 

1.012 

1.003 

0.977 

0.992 

-3.15 

0.990 

0.995 

1.001 

1.001 

0.986 

0.994 

-2.85 

0.977 

0.981 

0.987 

0.998 

1.006 

1.001 

-2.55 

0.981 

0.982 

0.978 

0.991 

1.031 

1.008 

-2.25 

1.003 

1.009 

0.982 

0.988 

1.052 

1.022 

-1.95 

1.030 

1.032 

1.002 

0.993 

1.058 

1.022 

-1.65 

1.048 

1.035 

1.025 

1.008 

1.037 

0.997 

-1.35 

1.034 

1.013 

1.029 

1.017 

0.982 

0.953 

-1.05 

0.971 

0.949 

0.986 

0.988 

0.893 

0.869 

-0.75 

0.852 

0.841 

0.880 

0.901 

0.772 

0.757 

-0.45 

0.684 

0.700 

0.721 

0.746 

0.631 

0.631 

-0.15 

0.496 

0.523 

0.536 

0.553 

0.487 

0.499 

0.15 

0.332 

0.339 

0.357 

0.367 

0.356 

0.370 

0.45 

0.197 

0.204 

0.218 

0.211 

0.244 

0.268 

Table  VI  (Cont.) 


Si(6;2)  Si ( 2 ; 2 ) Si(2;l) 

z 


n (z) 
e 

"hU) 

ng(z) 

nh<z) 

n (z) 
e 

nh(2) 

0.75 

0.103 

0.106 

0.122 

0.109 

0.158 

0.181 

1.05 

0.053 

0.055 

0.063 

0.050 

0.095 

0.122 

1.35 

0.026 

0.028 

0.031 

0.021 

0.057 

0.073 

1.65 

0.008 

0.010 

0.015 

0.008 

0.036 

0.048 

1.95 

0.002 

0.003 

0.006 

0.003 

0.019 

0.030 

2.25 

0.001 

0.001 

0.002 

0.001 

0.008 

0.016 

2.55 

0.000 

0.000 

0.001 

0.000 

0.003 

0.008 

2.85 

0.000 

0.000 

0.000 

0.000 

0.001 

0.004 

Table  VII  - Results  from  self-consistent  calculation  for  the  surface  properties 

of  EHD  in  Ge(4;2),  Ge(l;2),  Ge(l;l),  Si(6;2),  Si(2;2),  and  Si(2;l).  a is  the 

surface  tension  of  EHD.  p and  p,  are  the  bulk  chemical  potentials  of  elec- 

e h 

trons  and  holes.  These  are  obtained  from  the  fully  self-consistent  calculation 
(26 ) 

of  Vashishta  et.  al.  A$  denotes  the  dipole  layer. 


SYSTEM 

a (erg/cm  ) 

p (meV) 
e 

Uh(meV) 

V5e 

2A*(meV) 

SIGN  OF  CHARGE 

Ge(4;2) 

3.7xl0-4 

-3.62 

-2.21 

1.41 

0.77 

NEGATIVE 

Ge(l;2) 

l.OxlO-4 

-1.75 

-2.76 

-1.01 

-0.98 

POSITIVE 

Ge(l;l) 

0.2xl0'4 

-2.17 

-0.93 

1.24 

1.20 

NEGATIVE 

Si(6;2) 

87.4xl0"4 

-14.01 

-7.84 

6.17 

3.34 

NEGATIVE 

Si( 2 ;2 ) 

32.8xl0-4 

-8.87 

-9.89 

-1.02 

-1.02 

* 

NEUTRAL 

Si ( 2 ;1 ) 

11.4xl0-4 

-9.77 

-4.88 

4.89 

3.86 

NEGATIVE 

In  Si(2;2),  we  find  that  the  difference  in  the  chemical  potentials  of  holes 
and  electrons  ( = w - u - 2A<>)  is  a very  small  negative  quantity,  which  im- 
plies that  the  drop  is  positively  charged.  However,  considering  the  limit 
of  accuracy  of  the  calculation  we  can  only  conclude  that  the  EHD  is  neutral 
in  Si(2:2). 


FIGURE  CAPTIONS 


Fig.  1; 

Self-consistent  density  profiles  of  electrons  and  holes  in  Ge(4;2). 

Solid  curve-electron  density,  ng(z);  dashed  curve-hole  density,  n^(z).  The 
distance  z is  measured  in  units  of  excitonic  Bohr  radius,  a^C2  -77  X in  Ge), 

_ 3 

and  the  densities  are  in  units  of  a^  . In  Fig.  (1),  ng(z)  and  n^(z)  are 

3 

normalized  to  the  mean  particle  density  ( = 3/4Trr  where  r is  the  average 

s s 

interparticle  separation)  in  homogeneous  EHL.  Note,  the  hole  density  spills 

out  more  than  the  electron  density,  because  the  binding  energy  of  holes  is 

less  than  that  of  electrons  (y  - y =1.41  meV).  This  is  why  the  EHD  sus- 

h e 

tains  a negative  charge  in  Ge(4;2). 

Fig.  2: 

Normalized  self-consistent  electron  and  hole  densities  as  a function  of 

z in  Ge(l;2).  Solid  curve-electron  density;  dashed  curve-hole  density.  In 

Ge(l;2),  the  binding  energy  of  electrons  is  smaller  than  that  of  holes  (u^- 

y = -1.01  meV),  causing  the  electron  density  to  spill  out  more  than  the  hole 
e 

density. 

Fig.  3: 

Normalized  self-consistent  electron  and  hole  density  profiles  in  Ge(l;l). 
Solid  curve-electron  density;  dashed  curve-hole  density.  Note,  the  hole  den- 
sity tail  is  longer  than  that  of  electrons.  It  is  a consequence  of  the  fact 
that  the  binding  energy  of  electrons  is  larger  than  that  of  holes  (ji^  - y « 
1.24  neV). 

Fig.  4; 

Self-consistent  electron  and  hole  densities  as  a function  of  z in  Si(6;2). 
Unit  of  distance  is  excitonic  Bohr  radius,  ax  (=  49  8 in  Si),  and  the  unit  of 
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W F 


density  is  a, 


-3 


The  electron  and  hole  densities  are  normalized  to  the  mean 


particle  density  in  the  homogeneous  EHL.  Since  the  binding  energy  of  electrons 

is  larger  than  that  of  holes  (p,  - ii  = 6.17  meV),  the  electron  density  spills 

n e 

out  less  than  the  hole  density  in  Si(6;2). 

Fig.  5: 

Normalized  self-consistent  density  profiles  for  electrons  and  holes  in 
Si(2;2).  Solid  curve-electrons;  dashed  curve-holes.  In  Si(2;2),  the  binding 
energy  of  electrons  is  smaller  than  that  of  holes,  which  should  lead  to  a posi- 
tive charge  on  the  drop.  But  the  small  difference  in  the  chemical  potentials 
of  electrons  and  holes  (|jl^  - p | = 1.02  meV,  compared  to  6.17  meV  in  Si(6;2) 

and  4.89  meV  in  Si(2;l),  see  Table  VII)  is  almost  compensated  by  twice  the  di- 
pole layer.  Thus,  within  the  limits  of  accuracy  of  the  calculation,  our  con- 
clusion is  that  the  drop  is  almost  neutral. 


Fig.  6: 

Normalized  self-consistent  electron  (solid  curve)  and  hole  (dashed  curve) 
densities  as  a function  of  z in  Si(2;l).  It  is  evident  from  the  smaller  elec- 
tron density  tail  that  the  electrons  are  more  tightly  bound  than  the  holes, 

(y^  - 4e  = 4.89  meV),  thereby  giving  rise  to  a negative  charge  on  the  EHD. 


Figure  6 


